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Abstract 

This chapter is devoted to the recent theoretical results on the optical quantum control over 
charges confined in quantum dots under influence of phonons. We show that lattice relaxation 
processes lead to decoherence of the confined carrier states. The theoretical approach leading 
to a uniform, compact description of the phonon impact on carrier dynamics, perturbative in 
phonon couplings but applicable to arbitrary unperturbed evolution, is described in detail. Next, 
some applications are presented: phonon damping of Rabi oscillations in quantum dots and 
phonon-induced error of a single-qubit gate for an excitonic quantum dot qubit as well as for a 
semiconductor quantum dot spin qubit operated via a STIRAP transfer. 

1 Introduction 

With the state-of-the-art experimental techniques it is possible to control the quantum state of charge 
carriers in a quantum dot (QD). Many effects known from quantum optics of natural atoms have 
been demonstrated in these man-made systems: controlled coherent dynamics in these structures 
has been induced [1], Rabi oscillations have been observed [2-6]) and entanglement between states 
of interacting dots [7] has been demonstrated. 

However, unlike their natural counterparts, the artificial atoms are solid-state structures, embed- 
ded in the surrounding macroscopic crystal. Therefore, even in top-quality samples, some perturbing 
interaction effects are inevitable. The mutual influence of lattice deformation (phonons) and charge 
distributions is one of these inherent effects. There are three major mechanisms of carrier-phonon 
interaction [8]: (1) Coulomb interaction with the lattice polarization induced by the relative shift of 
the positive and negative sub-lattices of the polar compound, described upon quantization by longi- 
tudinal optical (LO) phonons; (2) deformation potential coupling describing the band shifts due to 
lattice deformation, i.e. mainly longitudinal acoustical (LA) phonons; (3) Coulomb interaction with 
piezoelectric field generated by crystal deformation (LA and transversal acoustical, TA, phonons). 
The lattermost effect is weak for globally charge-neutral excitations (e.g. excitons) in some systems, 
like InAs/GaAs, but may be of more importance for uncompensated charge distributions (e.g. ex- 
cess electrons) or for the properties e.g. of GaN dots [9,10], where charges are separated by large 
built-in fields. 

The coupling to the lattice degrees of freedom manifests itself in many ways in the spectroscopic 
properties of QDs. Resonant interaction with LO phonons leads to the extremely pronounced spec- 
trum reconstruction (formation of resonant polarons) [11-14], acoustic and optical phonons provide 
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a relaxation channel for carriers [15-17] with an important role played by phonon anharmonicity 
[18-21], phonon replicas and phonon-assisted transitions are magnified due to resonant interaction 
with quantized carrier states [22-27]. 

The perturbing effect of lattice modes has also been observed experimentally as a fast (~ 1 
ps) partial decay of coherent optical polarization induced by an ultra-fast laser pulse [28,5]. The 
theoretical analysis both in the linear regime [29,30] and in the nonlinear case [31,32] shows that this 
decay should be viewed as a trace of coherent lattice dynamics (due to lattice inertia) rather than as 
an effect of typical noise. The agreement between the theoretical modeling of carrier-phonon kinetics 
[32] and the experimental results [28] shows that phonon-related effects play the dominant role in 
the kinetics of a confined system on picosecond timescales. Phonon effects are also likely to play the 
leading role in the damping of Rabi oscillations induced optically in a QD [33-35] (although in the 
current experiments other effects are also important). Due to strong reservoir memory on timescales 
relevant for these processes, they cannot be fully understood within the Markovian approximations. 
In particular, the idea of a "decoherence time", with which the control dynamics competes, is 
misleading [36]. 

Apart from the general scientific interest, a strong motivation for studying the phonon processes 
taking place in semiconductor nanostructures comes from the recent proposals for defining a quantum 
bit (qubit, the basic unit of quantum information) in terms of orbital (charge) degrees of freedom 
of an exciton confined in a quantum dot [37] or in terms of spin states controlled via conversion to 
orbital degrees of freedom [38-43] . The recent experimental demonstration of a quantum logic gate 
operation with charge degrees of freedom in a QD [44] has strengthened this motivation even more. 
Understanding phonon-induced decoherence processes is essential for practical implementation of 
these novel ideas and currently seems to be one of the most topical issues in the field. 

In order to study the phonon effects for the general coherent-control and quantum information 
processing schemes, a theoretical method going beyond the perturbative treatment of the external 
driving is needed. In this chapter, we present one of such methods, treating the Coulomb interactions 
between the confined carriers and the coupling to the driving field exactly (non-perturbatively), 
while the phonon coupling is included as a perturbation. Although the range of validity of such a 
strictly perturbative treatment may be narrower than that of more sophisticated methods (e.g. the 
cumulant expansion technique [33]), the advantage of the proposed approach is that it yields closed 
formulas leading to a clear physical interpretation of the phonon decoherence effects. 

The chapter is organized as follows: The next section describes the system under study, intro- 
duces its model and discusses the unperturbed evolution under special conditions. In the section 
3 we describe the various carrier-phonon coupling mechanisms for single carriers and for excitons. 
The general idea of decoherence due to lattice relaxation is introduced in the section 4. Next, in 
the section 5 we develop the theoretical treatment for analyzing the phonon-related decoherence for 
an arbitrary system evolution. This is then applied to the description of phonon damping of Rabi 
oscillations (section 6) , optimization of control for an excitonic QD qubit (section 7) and fidelity of 
a spin qubit operated via STIRAP transition (section 8). The final section concludes the chapter. 

2 The system and the model 

The Hamiltonian of the system. We will consider a system of charges confined in a QD, 
coupled to coherent electromagnetic field and interacting with phonons of the surrounding crystal 
medium. The Hamiltonian for the system is 

H = H C + H ph + H int . (1) 

The first term Hq describes the carrier subsystem together with the external driving field. Through- 
out this chapter, we will assume that this driving field is strong enough and its quantum fluctuations 
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may be neglected, so that the field may be modeled classically. Usually, for a confined system only 
the lowest part of the discrete spectrum is relevant for reasonable pulse durations and intensities, 
so it is often convenient to assume that the total wavefunctions of the interacting carrier system 
are known (e.g. from numerical diagonalization [30]) and to express the Hamiltonian in the basis of 
these eigenstates. Thus, unless explicitly stated otherwise, we will denote the crystal vacuum by |0) 
and the confined exciton states by |n), n > 1. 

Numerical calculations for the harmonic confinement model [30] show that the lowest exciton 
states may be well approximated by products of the electron ground-state harmonic oscillator wave- 
function and a hole wavefunction corresponding to one of the low harmonic oscillator levels. This is 
due to the large hole effective mass compared to that of the electron. For the same reason, even if 
the geometrical confinement is the same for both particles, the Coulomb interaction shrinks the hole 
wavefunction while the electron wavefunction is only slightly modified (the electron excitation energy 
is much larger than Coulomb energy in a typical self-assembled structure). The separation between 
the ground state and the lowest dark excited state (hole excitation) is usually of several meV, while 
the closest bright state (the lowest one with the electron in an excited state) lies approximately at 
70 meV. 

Also, for the structure of excitonic levels as discussed above, the single-exciton assumption is 
reasonable under excitation with polarized light. Even though excited single-exciton states may 
be close to the ground state, they are formed by exciting the hole (which is much heavier than 
the electron), while the electron remains roughly in its ground single-particle state. Creation of a 
bi-exciton requires much higher energy, sufficient to promote both carriers to excited states. 

The second term in (1) is 

H ph = Y,^b{b k (2) 
fc 

and describes the energies of phonons, b^b^ being the phonon creation and annihilation operators 
(branch index will always be included into k, unless explicitly written). Since only long- wavelength 
phonons are effectively coupled to carriers confined in a QD (see discussion in the next section) we 
will here always assume linear and isotropic dispersion for acoustic phonons, so that uj s (k) = c s k, 
where c s is the speed of sound for the branch s (s = 1 for LA, s = t for TA). Also, because the 
dispersion of the LO phonons is weak around the zone center, the LO phonons (s = o) will be 
assumed dispersionless, uj q = Q. We will neglect any anharmonicity effects and assume the free 
harmonic evolution of the lattice subsystem in absence of the carrier-phonon coupling. Only the 
bulk phonon modes will be included. Although a quantum dot implies certain electrostatic and 
mechanical discontinuity of the system and the role of confined phonons may be discussed [18,45], 
invoking bulk phonons is usually sufficient for explanation of experimental features, like e.g. polaron 
resonances [11,12]. 
The last term, 

flint = E \n)(n'\F nn ,(k)(bi + b- k ), F nn ,(k) = F* n ,(-fc), (3) 

fe;n,n'>l 

is the exciton-phonon interaction expressed in terms of the discrete confined carrier states \n). 
Although inter-band phonon terms may appear e.g. as a consequence of strain-dependent inter- 
band couplings, in the following sections we deal only with intraband phonon effects, so that the 
interaction Hamiltonian conserves the number of quasiparticles of each kind. 

The phonon wavevectors are restricted to the first Brillouin zone. However, modulations of the 
band structure with spatial periods much smaller than the size of the carrier wavefunction cannot 
influence the carrier energy. In the harmonic oscillator approximation for carrier confinement, this 
leads to exponential cut-off of carrier-phonon interaction (explicit formulas are given in the section 
3 below) at the frequency ujq ~ c/l, where c is the sound speed, I is the carrier confinement size. 
Thus, the exciton is effectively coupled only to the long-wavelength part of the phonon modes. On 



3 



the other hand, vanishing of the coupling in the limit of k = reflects the insensitivity of the 
system properties to shifting the lattice as a whole. For gapless bosons with linear dispersion, the 
characteristic frequency determines the reservoir memory times and sets up the timescale of non- 
Markovian effects, related to "dressing" of the localized carriers with coherent lattice deformation 
field (see section 4). 



Unperturbed resonantly driven evolution Assuming that the laser field couples resonantly 
only to one transition (say, from the "empty dot" state |0) to the ground exciton state |1)) the 
Hamiltonian for the carrier subsystem may be written 

H C0 = 5»Xn| + \f(t) (e-*|l)(0| + e--*|0)(l|) (4) 

n 

(the energies E n are defined with respect to the "empty dot" state). 

Using the canonical transformation to the "rotating frame" , defined by the unitary operator 



exp 



— iuj \n){n\t 

n>l 



, U vw \n) = e' luJt \n), n > 1, 



we obtain from (4) 

H C = UlH C0 U rw - ]T \n)(n\ = - ]T A n \n)(n\ + \f{t) (|1)(0| + |0)(1| 



(5) 



n>l 



n>l 



where A n = uj — E n . In the special case of strictly resonant coupling to the ground-state excitonic 
transition, to = E±, so that the detuning from the ground-state transition vanishes, Ai = 0, and 
— A n = e n where e n are the intra-band excitation energies for a single exciton. 

Let us denote the evolution operator for the exciton-light system (generated by He, without 
phonon interactions) by Uo(t,s), where s,t are the initial and final times, respectively. In the 
special case of resonant coupling only between |0) and |1) state the evolution operator may be found 
explicitly (note that in this resonant case He commutes with itself at different times) . The result is 



U c (t,a) = oo8^(|0X0| + |1)(1|) - ,sin^(|0)(l| + |1)(0|) + ^ \n)(n\e— «~ s \ 



(6) 



n>l 



where 



m = f dr/(r) 

J s 

is the rotation angle on the Bloch sphere up to time t. 



3 Interactions of confined carriers with phonons 

For completeness, in this section we summarize the derivation of the coupling constants between the 
bulk phonon modes and the confined carriers in a polar and piezoelectric semiconductor [8,46,47]. 

The deformation potential. Any crystal deformation leads to shifts of the conduction (c) and 
valence (v) bands which are, to the leading order, proportional to the relative volume change. The 
corresponding contribution to the energy of electrons (e) and holes (h) in the long-wavelength limit 
is, therefore, 

ff (DP) _ , AP . _ Sv 
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where a e /^ are the deformation potential constants for electrons and holes and V is the unit cell 
volume. Using the strain tensor a, 



1 / dui duj 
Tij = 2 I c>~ + 9r7 



one may write 

H ^/h ] = T^ c /hTra = Tff c/h V • u(r), 
where u(r) is the local displacement field. The displacement is quantized in terms of phonons, 



»M=^E^2^^(^ + ^)e*', (7) 



where uj(k) is the frequency for the wavevector k, e^^ = — e_fc s is the corresponding real unit 
polarization vector, and p is the crystal density. Only the longitudinal branch contributes to V ■ u 
in (7) and the final interaction Hamiltonian in the coordinate representation for carriers is 



<r^,^£/^K.+^,)^. (b> 

In the second quantization representation with respect to the carrier states this reads 
H iT = EH^WKKW = E I™ d\r n {r)V{r)AAr)aia n , 

, , J — oo 

nn' nn' 

4 E ^n'f^Uk) + bl^) . (9) 



^ fan' 



where 



with the formfactor 



C (DP) . / hk 



/oo 
d 3 rr n (r)e^ n ,(r) = ^ n (-fe). (10) 
-oo 

The interaction Hamiltonian (9) conforms with the general form assumed in (3). 

General properties of the formfactors. While the common coefficient of the coupling Hamil- 
tonian contains the fundamental and material-dependent constants and reflects the general electrical 
and mechanical properties of the semiconductor system, the formfactor (10) contains the informa- 
tion about the geometry of the confinement and the resulting properties of wavefunctions. In this 
sence, it is the "engineerable" part of the carrier-phonon coupling. 

Prom orthogonality of single-particle states one has immediately J- nn >(0) = S nn i. If the wave- 
functions are localized at a length I, then the extent of the formfactor is ~ l/l. Thus, for carrier 
states localized in a QD over many lattice sites and smooth within this range, the functions T nn '{k) 
will be localized in the fc-space very close to the center of the Brillouin zone. 

As an example, let us consider the ground state of the harmonic oscillator potential, 



1 fr±\ 2 1 [ z^ 2 



2\h 2\l 



(11) 



where r± is the position component in the xy plane and l±,l z are the localization widths in-plane 
and in the growth (z) direction. The corresponding formfactor is then easily found to be 



J-'(k) = exp 



k±l±\ 2 1 (k z l 



2' 

'°z b z 



(12) 
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Piezoelectric interaction. A propagating phonon wave in a piezoelectric medium induces a 
polarization field which affects the carriers by means of the Coulomb interaction. If the crystal 
deformation is described by the strain tensor a then the piezoelectric polarization is P = da, where 
d is the piezoelectric tensor. 

From the Maxwell equation for plane- wave fields (in absence of external charges and currents), 

ik ■ E = ——ik ■ P, ik x E = iuB, 



ik B = 0, 



c 2 ik x B = -iujE - —itoP, 



one has 



c 2 -[k(k ■ E) - k 2 E] = -uj[E + —e kE + — da], 

LO ' Eq Eq 



where we have used the relation P = eqkE + da. 
For the transversal component we get 



c 2 -k 2 E A 

LO 



-LO 



1 



E ± (1 + k) + — (da) ± 



hence, 



For strain fields associated with phonon propagation, to and k are phonon frequency and wavevector 
and Lo/k = c s <C c so that the transversal component vanishes. Thus, the piezoelectric field ac- 
companying a phonon wave is purely longitudinal and one may introduce the appropriate potential. 
One has explicitly for the longitudinal component 



±eEu = T 



(da) || 



where 



V(r] 



Eo 1 + K 

±e (do-) || 
= -i- ^ 

K E§E\ 



W(r), 



Using (7), one finds the strain tensor 

1 1 



E 



-(b s ,k + bf,-k) [(cs,k)ikj + {e a , k )jki] e 



ikr 



2 ^ iX V 2M ^s, k 

In the zincblende structure, one has 

dxyz = dy ZX = d ZX y = d, dijk = dikj , 
and the other components vanish. Hence 

(<i<7)|| = 2—(k x a yz + k y a zx + k z a xy ) 



and 



H 



(PE) 



=R 



n 



de 



M s (k)(b kyS + b>)e 



„ik-r 



where the polarization-dependent geometrical factor is 

M s (k) = 2 k x ky{e s ^) z + k y k z (e s ^) x + k z kx(e s ,k)y 



(13) 
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For the choice of the phonon polarizations (1-longitudinal, tl,t2— transversal) 



e\,k = = (cos0cos0, cos0sin0, sin0), 
e t i,fe = (- sin 0, cos 0,0), 
et2,fc = (sin 6 cos 0, sin 9 sin 0, — cos 0) , 



the functions M s are 



Mi(0,0) = ^ sin 20 cos sin 20, (14) 

M u (0,0) = sin 20 cos 20, (15) 
M t2 (0,0) = (3 sin 2 - 1) cos sin 20. (16) 

Finally, in the second quantization representation, 

< h E) = E^r^wn^M (17) 

= 7^ E 4«n'/ e ( / h E L', s ( fc ) ( fo fc, S + & -fcJ > 
* k,s,nn' 

where 

= ^Vra)^ M " ( ^""' (fc) ' 

with the formfactor given by (10). 

Note that the LA phonons are coupled to carriers both by the deformation potential and by the 
piezoelectric coupling, while the TA phonons only by the latter. The piezoelectric coupling results 
directly from the Coulomb interaction and has exactly opposite value for the electron and the hole. 
In contrast, there is in general no fixed relation between the values of the deformation potential 
constants. 



Frohlich coupling An LO phonon propagating in a polar medium is accompanied by a polar- 
ization field resulting from the relative shifts of the positive and negative ions forming the crystal 
lattice. Similarly as in the case of piezoelectric coupling, one shows that the resulting electric field 
is longitudinal and may be associated with a potential. Both the polarization and electric field, 
and hence the potential, are proportional to the LO phonon displacement. The derivation of the 
proportionality constant may be found in [8,46]. The resulting interaction energy for the charge 
carrier at point r is 

' x-» t I liV. 

k 

where e = (1/eoo — l/eg)" 1 is the effective dielectric constant. 
In the occupation number representation this reads 



H^ = T ^y:^ ^L(b k + bl (is) 



e/h 



L E 4«n'/ e ( / P h E L(fc) {b k ,i + blj) , (19) 



^ knn 



where 



{F[) _ e i_hV_ 

/e/h,nn' -=F fc y 2v e e [ 
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Phonon coupling for exciton states Above, we have derived the interaction Hamiltonian in 
the single-particle basis. However, most of the following deals with excitonic states, i.e. states of 
confined electron-hole pairs interacting by Coulomb potentials. 

Both carriers forming the exciton couple to phonons according to (9), (17) and (19). If the lowest 
exciton states are assumed to be approximately of product form, as discussed above, i.e. 

\ n ) = 4,i a h,J°>> 

then one gets from (9) the following coupling constants for the deformation potential interaction in 
the excitonic basis 

FS V \k) = ] j^ i {^(k)-a h ^(k)), (20) 

Fif P \k) = -]f^^\k), 1>1. (21) 

It is essential to note that, due to different deformation potential constants a e /^, none of these 
couplings vanishes even if the electron and hole wavefunctions are the same, leading to identical 
single-particle formfactors. 

For the Frohlich coupling to LO phonons one has 

= (>*?<*> < 22 > 

The form of the coupling constants for the piezoelectric interaction is analogous. In the case of 
Frohlich and piezoelectric coupling, the "diagonal" coupling Fu vanishes if the wavefunctions overlap 
exactly. This results from the fact that both these interaction mechanisms are related directly to the 
Coulomb interaction between the confined charge distribution and the phonon-related polarization 
field. For carriers localized in the same spatial volume but overlapping only partly there will be still 
some cancelation effect. It is remarkable, however, that even for perfect cancelation of the diagonal 
coupling, the other contributions F± n , n > 1 are not decreased. As we will see, this will lead to 
additional effects for LO phonons, while in the case of piezoelectric coupling these contributions are 
negligible due to large energy level spacing compared to the acoustic phonon frequencies. 
The spectral properties of the lattice are characterized by the phonon spectral densities 

Rnn>, mm >(u) = ^IneH + E F nn'(k)F^>(k) [6(<j - u s (k)) + 8{u + u s (k))\ , (24) 

k,s 

where hb(w) = — ne(— oj) — 1 is the Bose distribution function. These functions depend on the ma- 
terial parameters and system geometry and fully characterize the properties of the lattice subsystem 
at the level of perturbation treatment discussed in this chapter. As will be seen later, they are one 
of the two "building blocks" for the perturbative calculation of the phonon effects on a quantum 
evolution. Here, we will need only a subset of these functions, Ri(uj) = Rii,ii(us). 

Let us consider the functions i?| LA ^ (lo), corresponding to the LA phonons. Here the total coupling 
F^ LA ^ contains both DP and PE contributions. However, under reasonable symmetry assumptions 
the formfactors may be chosen either real or purely imaginary [see e.g. the explicit form for the 
harmonic confinement, Eq. (12)], while the coefficients in (9) and (17) are purely real and purely 
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imaginary, respectively. Thus, the spectral densities split into two independent contributions. Sub- 
stituting (20,21) into (24) and performing the summation over k in the usual continuum limit 

\ ( 2 -) 2 J 

one obtains 

rI DP \lo) = R DP io 3 [n B (io) + l]/,( w ), (25) 

where 

4n 2 hpcf 

and fi(cu) are certain functions that depend on the wavefunction geometry, having the property 
/i(0) = 1 and fi(0) = 0, I > 1. A similar procedure leads to an analogous result for the piezoelectric 
contribution from LA phonons. 

In the case of dispersionless LO phonons 0J o (k) = O and the spectral densities are 

R^ t \uj) = Rptlnsiu) + l\[S{u - n) + S(u + Q)] 
R^ito) = R^mlnBicj) + l\[6(u - Q) + 8(co + n)] 



where 



r 30e 2 {ll - _ tte 2 

-"Fr ~ n 7^= — ^7 TcTi ' -"Fr 



SnV^eolL 16L 4 ' tr 87rV27re eZh' 

and 77(00) = 1) ^(oi) = 1/4, ??(02) = 3/32, 7710 = 3/16. Here we numbered the states by (N, M), where 
M is the total angular momentum and N is another quantum number, and introduced the averaged 
carrier localization L 2 = (l 2 + l 2 J/2, where l e ,l\i correspond to l± in Eq. (11) for electrons and holes, 
respectively, and / 2 — l\ <C I? . 



4 Decoherence of carrier states by a dressing process 

Any optical experiment or optical control of carrier states in a quantum dot is based on the coupling 
between the carriers and the electric field of the electromagnetic wave. By using high-power laser 
pulses this coupling can be made strong, leading to fast excitation of confined excitons. However, 
this coupling involves pure electronic degrees of freedom, inducing only carrier transitions which 
must then be followed by lattice relaxation to the potential minimum corresponding to the newly 
created charge state. In the case of the optical exciton creation, the initial lattice configuration 
corresponds to zero polarization field and zero deformation, while the eigenstate of the interacting 
system involves some lattice polarization (LO phonon field) and some deformation (LA phonon 
field). Thus, the process of the lattice relaxation may be viewed as the formation of the phonon 
dressing around the confined exciton. The lattice relaxation process, accompanied by a non-zero 
probability of phonon emission, only takes place if the result of a control pulse is a one-exciton state. 
Hence, if the pulse creates a superposition, the lattice dynamics is different for the two components 
of the superposition state, leading to entanglement between the carrier and lattice degrees of freedom 
and to decoherence. Physically, detection of a phonon carrying out the excess energy during the 
lattice relaxation process reduces the superposition of carrier states to the one-exciton state, since 
the other state (empty dot) is never accompanied by the lattice relaxation. It should be kept in 
mind, however, that unlike in the phonon-assisted carrier relaxation processes, the average number 
of phonons emitted during the dressing is much less then one (in the weak coupling case) and even 
in the limit of t — > 00 the emission probability is below 1, leading only to partial loss of the carrier 
coherence [29,30]. 



9 



In this section we discuss the carrier-phonon kinetics after an ultra-fast excitation from the point 
of view of decoherence and "information leakage" from the carrier subsystem [30] (see [29,31,32] for 
an exhaustive discussion from the point of view of the optical properties). 

We consider a quantum dot in which an exciton is created by a very short laser pulse. We assume 
that the pulse is much shorter than the periods of the phonons but long enough to ensure spectral 
overlap only with one (ground-state) exciton transition. We allow, however, for the existence of 
dark states of energy much lower than the closest bright state, as discussed in the section 2. In 
accordance with the discussion in section 3, we include in the interaction Hamiltonian (3) the 
deformation potential coupling to LA phonons and Frohlich coupling to LO phonons. 

In order to quantify the quality of a manipulation on the quantum state from the point of view 
of the desired goal, one defines the fidelity as the overlap between the desired final state obtained 
from the initial state |-0o) by the unperturbed evolution Uq and the actual state described by the 
reduced density matrix of the carrier subsystem p(t), 

F = 1 - 8 = (Muh(t)pWc(tMo)T, (26) 

where (-)t denotes thermal average over the initial lattice states. If the initial carrier state is the 
vacuum |0) and the operation, performed at t = 0, is an ultra-fast rotation by the angle a on the 
Bloch sphere corresponding to the |0) and |1) states than the unperturbed state for t > is 

\m) = U c (t)\0) = cos||0) -i e -^sin||l). (27) 

Let us calculate the actual state, including phonon perturbations, in this ultrafast limit. The 
diagonal elements of the density matrix remain constant, since phonon processes cannot change the 
excitonic occupations, 

Poo = cos -, pu = sin -, 
while for the non-diagonal ones one finds 

pio(t) = Tr [f?(t)|0Xl|] = ~ismaG(t), (28) 

where we have substituted p(0 + ) = |^(0))(^(0)| <g> u>l, with | -0(0)) given by (27), and defined the 
correlation function 

G(t) = (u\t)\0)(l\U(t)\l)(0\) T , 

where U(t) is the evolution operator for the total interacting system. Using the definition (26), the 
fidelity loss 5 may be expressed as 



S(t) = - sin 2 a [l - Re (G{t)e iElt )] . (29) 



2 

In order to find the correlation function we describe the carrier-phonon kinetics after an abrupt 
excitation using the basis of perturbative eigenstates of the Hamiltonian H = Hq + H irA [20,14]. 
Thus, we introduce the new states and operators 

\n) = e s \n), (3 k = e s 6 k e~ s , 

corresponding to dressed particles. The system Hamiltonian expressed in terms of these operators 
is diagonal up to second-order corrections if the anti-hermitian operator S is chosen as 

S= E * ( :l(k)\n)(n'\(b k , s -bl ktS ), 

n,n'>l;k,s 



10 



where 



E n / -£ n + ftu; s (fc)' 



(s) 

and F„ n , are given by (20-23). Therefore, the evolution of these new states and operators is free, 



a n {t) = \n) = e-^ln), ft, fc (t) = ^ k e~^\ 

where E n is the exciton energy including corrections resulting from the interaction. 

To the lowest order, the creation operator for the ground exciton state may now be written as 



|1)(0| = e S |l)(0| e - s 

= !"^E \<t>&(k)\ 2 Wl, k Ps, k + 1) I |1)(0| 



n.s.k 



Using these formulae, it is easy to write down the correlation function 
G(t) = (l-l5:i45(fc)| 2 (2n s , fc +l)| e -^ 



N 

n,s,k 

4E^+*2wi 2e " (l " /SK( ' !),i . 



N , 

n,s,k 

where n s ^ are the thermal occupation numbers for phonon modes. It should be noted that the 
perturbative treatment is valid only when the above correlation function remains close to one. For 
a fixed strength of interaction this requires a low enough temperature (in practice, for a typical 
InAs/GaAs dot, this means T < 100 K). 

In order to explain the time evolution of this function, let us remember that the functions $^,(fe) 
contain the formfactors (10), effectively selecting a certain wavenumber range k ( Q s) ±±Afc( s ), centered 
around k^ for each branch of phonons (Ak^ ~ l/l, where I is the dot size). This corresponds to 
a certain frequency range ui^ ± ^Aw' s ', around a central frequency lu^ . The phases that enter in 
the summation in (30) at the time t spread effectively over the angle range Uq i} t ± iAw (s) t. When 

_2tt 

asymptotic value of the corresponding error is 



this phase spreading reaches 2tt, i.e. for t ~ . 2 T s) , the last two terms in (30) become small. The 



^ = ^ sin2 «^Ei4S(fc)i 2 (2^ fc + i). 

k,n,s 

In the case of acoustical phonons this value critically depends on temperature (Fig. la). The 
decoherence time depends on the phonon dispersion (the frequency range Auj^) and weakly evolves 
with temperature. For the acoustical phonons, r ~ 2irl/c ~ 1 ps (Fig. la). For the nearly 
dispersionless optical phonons, the dynamics is dominated by slowly vanishing coherent phonon 
beats (Fig. lb,c). However, in order to excite these oscillations one needs pulses shorter than LO 
phonon periods, i.e. of durations ~ 10 fs. Moreover, it is known that the anharmonic LO-TA 
interaction [48,49] (not included in the present description) acts on much shorter timescales and 
may be expected to considerably shorten this time. 
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t [ps] t [ps] t [ps] 

Figure 1: Dressing- induced decoherence of the exciton due to LA phonons (a) and due to LO 
phonons: oscillations with LO phonon frequency (b) decay on a very long timescale (c) (anharmonic- 
ity effects are not included). Four lowest exciton states calculated by numerical diagonalization [30] 
were taken into account. 




Figure 2: Time evolution of the mean lattice deformation (acoustic phonons) at the dot axis (x = 
y = 0). The time-independent deformation, corresponding to the coherent phonon dressing, is 
formed around z = (i.e. in the dot area) within ~ 1 ps. This is accompanied by emission of a 
phonon packet carrying away the excess energy at the speed of sound, seen on the plot as the ridge 
and the valley with growing distance from the dot. 



It is possible to see directly that the described process consists in the formation of a coherent (i.e. 
with non-vanishing mean displacement) phonon field corresponding to the classical lattice deforma- 
tion around the confined charge. Let us consider the mean lattice displacement corresponding to the 
branch s at the point r after time t after rapid creation of an exciton. Using (7) and transforming 
to the dressed basis one obtains in the lowest order 

<iKM)|o> T = 1 E ^.JL-*^ [« w (i - *m*)] . 

After a sufficiently long time the oscillating term averages to zero (around r = 0) and a time- 
independent displacement field is formed. The Fig. 2 shows the mean lattice deformation due to the 
deformation potential coupling on the dot axis, r = (0, 0, z), at various instances of time. A simple 
and intuitive picture emerges: formation of the lattice deformation, corresponding to the displaced 
equilibrium, is accompanied by emitting a phonon packet that carries the excess energy away from 
the dot with the speed of sound. It is clear that as soon as these wavepackets leave the volume of 
the dot, some information about the carrier state is carried into the outside world. Also, since the 
only controlled interaction is that with the confined carriers, once the phonons have been radiated 
out the reversibility is lost. 

We have shown that the carrier-phonon coupling leads to entanglement of the carrier system and 
the surrounding lattice. This, in turn, reduces the degree of coherence of the carrier subsystem. This 
decoherence is related to the spontaneous relaxation of the lattice to the new equilibrium defined by 
the carrier-phonon interaction and results from the non-adiabaticity of carrier evolution with respect 
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to the lattice response times. It may be expected that for a slow enough evolution most of the lattice 
modes will follow adiabatically, thus creating the coherent dressing cloud in a reversible way. In the 
following sections we derive the formalism suitable for the description of the carrier-phonon kinetics 
and apply it to a few problems. 

5 Evolution of the density matrix for a driven carrier-phonon sys- 
tem 

In this section we derive the equations for the reduced density matrix of the carrier subsystem in 
the leading order in the phonon coupling, assuming that the unperturbed evolution is known. This 
somewhat lengthy and technical derivation yields very simple and intuitive formulas that may be 
easily applied to a range of problems to be discussed in the subsequent sections. 

We will consider a system composed of charge carriers (electrons and holes), localized in one or 
more quantum dots, driven externally and interacting with phonons, as described by the Hamilto- 
nian (1). As already mentioned, the evolution of the carrier subsystem is generated by the (time 
dependent) Hamiltonian He, describing the properties of the system itself as well as its coupling to 
driving fields. In this chapter, we deal only with optical driving (although not necessarily via direct 
excitation of dipole-allowed transitions), but other types of external fields may be treated on the 
same footing. The evolution of the phonon subsystem (reservoir) is described by the Hamiltonian 
H p h [Eq. (2)]. We will restrict the discussion to the free phonon evolution, neglecting anharmonicity 
effects. The evolution operator for the driven carrier subsystem and free phonon modes, without 
carrier-phonon interaction is 

U (t) = U c (t) ® e- iH ^- s \ 

where Uc(t) is the operator for unperturbed evolution of the carrier subsystem (we suppress the 
initial time in U(t,s)). 

The carrier-phonon coupling is described by the interaction Hamiltonian (3) which may be 
written in the form 

V = Y, S nn'®Rnn', (30) 
nn' 

where n, n' run over the carrier subsystem states, S nn / (possibly time-dependent) act in the Hilbert 
space of the carrier subsystem while the time-independent R nn > affect only the environment. It is 
convenient to allow non-Hermitian operators S nn > and R nn '] however, we demand the symmetry 
relation 

S nn ' = Sn'n-, R nn > = Rn'n (31) 

which guarantees hermicity of the Hamiltonian (30) and is also explicitly satisfied by (3). 
We will assume that at the initial time s the system is in the product state 

Q(s) = | VoX^o | ®PT, (32) 

where \ipo) is a certain state of the carrier subsystem and pt is the thermal equilibrium distribution 
of phonon modes. Physically, such an assumption is usually reasonable due to the existence of two 
distinct time scales: the long one for the carrier decoherence (e.g. 1 ns ground state exciton lifetime 
[28,50]) and the short one for the reservoir relaxation (1 ps dressing time [28-30]). 

The starting point is the evolution equation for the density matrix of the total system in the 
interaction picture with respect to the externally driven evolution Uo, in the second order Born 
approximation with respect to the carrier-phonon interaction [51] 

Qit) = ~Q(s) + ^[ dr[V(r), g(s)} - 1 J* dr £ dr'[V(r% [V(r'% g(s)}}, (33) 
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where 

g(t) = ui(t)g(t)U (t), V(t) = tf (t)VU (t) 

(it should be kept in mind that V may depend on time itself). 

The reduced density matrix of the carrier subsystem at time t is 

p(t) = U c (t)p(t)Ul(t), p(t) = Tr R g(t), 

where the trace is taken over the reservoir degrees of freedom. The first (zeroth order) term in (33) 
obviously gives rise to 

P (0) (t) = u c (tMo)(Huh(t) = IMWMQl- (34) 

The second term vanishes, since it contains the thermal average of an odd number of phonons. The 
third (second order) term describes the leading phonon correction to the dynamics of the carrier 
subsystem, 

~P {2 \t) = -± fjT jyr'TT R [V(r'),[V(r"),g(s)]]. (35) 
First of the four terms resulting from expanding the commutators in (35) is 

(i) = -WoXM 

where ^ 

Qt = IlXIXl / dT / dT ' S nn'{T)S mml {T r ){R nnl {T -T')R mml ). (36) 

ft i i J s J s 

nn mm 

The operators S and R are transformed into the interaction picture in the usual way 

Snn'(t) = tf (t)S nn ,U {t), Rnnit) = tf (t)R nn ,U (t) 

and (O) = Tt r [Opt] denotes the thermal average (obviously [Uo(t),pr] = 0). 
The second term is 

I rt rr 

( n ) = ^EE / dr / dT '\^o)(^o\S mm '(r')S nn/ {T){R mm/ (T' - T)R nn >). 

tl , i J s J s 

nn' mm' 

Using the symmetry relations (31) one has 

(Rmm'ir' -T)R nn >)* = {R n ' n {r - T r )R m i m ), [S mm > (r')S nn / (r)] * = S n i n (T)S m > m (T'), (37) 
hence this term may be written as 

(ii) = -woXmqI- 

In a similar manner, using the symmetries (31), the two other terms may be combined to 

(III) + (IV) = 6 t [IVoXV'ol] ■ 

where ^ 

®t \P\ = 72 E E / dT [ T dT'S nn/ {T') P S mm/ {T){R mm/ {T - T')R nn ,). (38) 

Tl , , J s J s 

nn' mm' 

In terms of the new Hermitian operators 

A t = Q t + Ql h t = ±- i (Q t -Q\), (39) 
the density matrix at the final time t (34,35) may be written as 

P (t) = u c (t) (iVoXV'ol - i [h, \MM - \ i A t, IMM + $tUo){M]) 4- ( 40 ) 
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The first term is a hamiltonian correction which does not lead to irreversible effects and, in principle, 
may be compensated for by an appropriate modification of the control Hamiltonian Hq. The other 
two terms describe processes of entangling the system with the reservoir, leading to the loss of 
coherence of the carrier state. 

Let us introduce the spectral density of the reservoir, 

Rnn'mm' M = ^2 / dt{R nn ,{t)Rmm>)e ibJt - (41) 

If the operators R nn ' are linear combinations of freely evolving bosonic modes (only this case will 
be considered in this chapter), 

Run' = Rl n = -!= E F nn>{k) (b h + bQ , (42) 

with F nn i{k) = F*, n (—k) (branch index implicit in fe), then (41) coincides with (24). 
With the help of (41) one may write 

nn' mm' 

where the frequency-dependent operators have been introduced, 

Y nn .(w) = fdrSnn^e™. (44) 

Using (41) one has also 

Qt = EE^ f dT r^(T-r')^(r)5 mm ,(r')i?n„',w^)e-^ {T - T,) . 

/ / J s J s 

nn mm' 

Next, representing the Heaviside function as 

r HtJ p-iu't 
9(t) = -e^ t J dU 



2m J - w + i0+ ' 



we write 



2m u' -uj + i0+ 



ft = -EE/^w/ i;Wl '-' M 

nn' mm' 

= _ EEy dujR n n',mm'(^>) 



nn' mm' 

dJ 



to 



where V denotes the principal value. 

Hence, the two Hermitian operators defined in (39) take the form 



nn' ,mm' 

{u)Y ] n , n {uj)Y mrn ,{u) (45) 



and 

dw'Yl ln {u')Y mml {u') 



h = E E / duRnn',mm'^)V J ^ Y ^ U f™'*>\ ^ 
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where we have used the relation 

R nn',mm'{u) = R m 'm,n'n{^), 



resulting from the definition (41) and the symmetry relation (37). 

Using the definition of the fidelity (26) and the Master equation (40), the error may be written 
in a general case as 



(V'ol^l^o) - (i>o $[|^oX^o|] 



Y>o>- 



It should be noted that the unitary correction generated by ht does not contribute to the error at 
this order. 

Using the definitions (38,45) this may be further transformed to 



5= J dujR rm , :mm ,(Lj) (ifo yJ, n (w)P- L y mm /(^) 4>o) , 

nn'mm' 

where P -1 is the projector on the orthogonal complement of |V>o) m the carrier space. 



(47) 



6 Phonon- induced damping of Rabi oscillations 

Introduction. Observation of Rabi oscillations of charge degrees of freedom confined in quantum 
dots [2-6] is believed to be a fundamental step towards quantum control of these systems. However, 
the coherent dynamics of the confined carrier states is very sensitive to any interaction with the 
macroscopic number of degrees of freedom of the outside world. In fact, so far it has always turned 
out that experimentally observable Rabi oscillations deviate from the ideal ones, the discrepancy 
being larger for stronger pulses. In principle, this might also be explained by experimental condi- 
tions or environmental perturbation: scattering by weakly localized excitons around an interface 
fluctuation QD (further confirmed by increasing decay for stronger pulses) [2], tunneling to leads 
in the photodiode structure (on ~ 10 ps timescale) [4], or dipole moment distribution in the QD 
ensemble [28]. 

One might believe that all the perturbation comes from sources that may be removed or mini- 
mized by technology improvement and by optimizing the experimental conditions and hence produce 
no fundamental obstacle to arbitrarily perfect quantum control over the excitonic states. However, 
in every case the QDs are inherently coupled to the surrounding crystal lattice. The recent theoreti- 
cal study [33] on optical Rabi flopping of excitons in QDs driven by finite-length optical pulses in the 
presence of the lattice reservoir shows that exponential damping models fail to correctly describe the 
system kinetics. The appropriate quantum-kinetic description yields much less damping, especially 
for long pulses. It turns out that the lowest "quality" of the Rabi oscillation is obtained for pulse 
durations of a few ps, while for longer durations the damping is again decreased. 

In this section we discuss both qualitative and quantitative explanation [34] of the mechanism 
leading to the phonon-induced damping reported in the theoretical and experimental studies. We 
show that the carrier-phonon interaction responsible for the damping of the oscillations has a res- 
onant character: While in the linear limit the system response depends only on the spectral de- 
composition of the pulse, the situation is different when a strong pulse induces an oscillating charge 
distribution in the system. In a semi-classical picture, this would act as a driving force for the lattice 
dynamics. If the induced carrier dynamics is much faster than phonon oscillations the lattice has 
no time to react until the optical excitation is done. The subsequent dynamics will lead to exciton 
dressing, accompanied by emission of phonon packets, and will partly destroy coherence of super- 
position states [28-31] but cannot change the exciton occupation number. In the opposite limit, the 
carrier dynamics is slow enough for the lattice to follow adiabatically. The optical excitation may 
then be stopped at any stage without any lattice relaxation incurred, hence with no coherence loss. 
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The intermediate case corresponds to modifying the charge distribution in the QD with frequencies 
resonant with the lattice modes which leads to increased interaction with phonons and to decrease 
of the carrier coherence (see Ref. [52] for a simple, single-mode model). 



The formalism In one version of the experiment [2,4,5] one measures the average occupation of 
the exciton ground state, (l|p(oo)|l), after a resonantly coupled pulse of fixed length but variable 
amplitude, starting with the system in the ground state |0). According to (27), in the ideal case the 
final occupation should vary as 



(l|p(oo)|l) idea i = sin 2 \a, 



where 



a = $(oo) 



drf(T), 



is the total pulse area. Within our formalism, the final occupation of the exciton state is given 
directly by the appropriate diagonal element of (40), once the spectral integrals entering in (38) and 
(45) for a given pulse are calculated. 

The interaction Hamiltonian has the form (3), with the coupling constants (20-23). Hence, 
Snn' = I^X^'I) n = 1, 2, . . . and, using the explicit form of the evolution operator (27), the operators 
(44) may be written in the form (n, n' > 1; we integrate by parts in order to extract the oscillatory 
contributions) 



Y-i 



li 



Y 



nl 



(loxi 
-(|i)(i| - |0)(0|) 

-|n)(0| 
-|n)(l 



smae 



I 

1 

ie n s 



cos ae iu!t -e iu!S -K^(u) 



(48) 



lo + e Tl 

i 



sin 



lo + e n 



cos r e i(«+en)t _ e i{u+e n )s _ K ^ 2 \lO + e r , 



y 



\n)(n'\ 



,-i{e„-e n ,)s 



lo + e n - e n / 

where a = $(oo) is the total rotation angle and 

rt ■ d 



J(u+t n -e n ,)t _ i{w+e n -e n ,)s 



kW>{w)= f dre tuJT ^-sm^{r), K^\oj) = f dre iulT ^- cos //$(r). 

J s J s 



(49) 



(50) 



(51) 



In the limit of s — > — oo, t — > oo (i.e. measurement after sufficiently long time compared to the 
reservoir memory), the functions (51) actually depend only on lot p , where r p is the pulse duration. 

Since the unperturbed evolution operator (27) has a block-diagonal structure and conserves 
the subspace spanned by |0) and |1), the value of pn(t) is determined by the matrix elements of 
pit) between these two states, while pu for I > 1 has the same form as in the interaction picture. 
Retaining only non-vanishing contributions, these matrix elements are, from (40), with |-0o) = |0), 



Poo = 1-E / dwRi(u)\(l\Y n \0)\ 2 , 
i>i 

Poi = J dujRxiuj) 



(52) 
(53) 



< ( (OlYnlOXllYnlO)* _ l(0|y n |0)*(0|yn|l) - ^l|Yii|0)*<l|Y n |i; 
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Figure 3: Phonon spectral density R(uj)/uj 2 (dotted) for deformation potential coupling to LA 
phonons at T = 10 K and the nonlinear pulse spectrum S a (uj) (solid lines) for r p = 1 ps, for rotation 
angles a as shown. 



l>2 



/ duR^HllYnWillYnll) 



Pu = J dc^HK^i|0)| 2 , I 



> 1, 



(54) 



where the short-hand notation Ri(lo) = Ruji(lo) has been introduced [the definition is given by 
Eq. (24)]. Apart from the dissipative correction described by the above formulas, there is an 
additional component, generated by the hermitian operator h t [cf. Eq. (40)], describing a unitary 
correction, e.g. light coupling renormalization and energy shifts, which may be canceled by an 
appropriate modification of Hq (these effects may lead e.g. to intensity dependence of the observed 
Rabi frequency [33]). 

The phonon spectral densities Ri{uj) are obtained from (24). We assume that no phonon-induced 
transitions are possible, i.e. the spectral densities i?;(±q) vanish. For any physical carrier-phonon 
coupling at arbitrary temperature one has also R(u) ~ uj n , n > 2 for uj — > 0. Thus, after substituting 
(48-50) into (52-54) and taking the limit s — ► — oo, t — > oo the contribution from the oscillating 
terms vanishes and the relevant elements of the density matrix, upon transforming back to the 
Schrodinger picture using (27), are 





2 a 


Poo 


= cos — 
2 




2 Q 


hi 


= sin — 
2 



did 



ili(w) 



E 

i>i 



Rl{u -ei) 
010 s o^[uj) 



US" 



+ / duo 



l>i J 



dio 



Rl{w — e l) 



l>i 



a 



where 



Pu = 
S c (u) 



sm 



cLo- 



R,.(uj 



+ / dw 



Ri{uj - ei) 



LO^ 



S c (oj) 

[s h (u) + s c (w) 



{cosa|^ s (1) (w)P 



sin aRe 
1 



cos 2 ^\KW 2 \u,)\ 2 - - sin aRe [K^l 2 > {uj)K^ 2 \u) 



sin 2 ^\K^ 2 \u)\ 2 + \ sinaRe [k^ 2 >(u;)K^ 



(55) 

(56) 
(57) 
(58) 

(59) 
(60) 
(61) 



The formulas (55-57) quantify the idea of resonance between the induced dynamics and lattice 
modes: in Fig. 3 the phonon spectral density is compared to the non-linear frequency characteristics 
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Figure 4: Pulse-area-dependent Rabi oscillations for various pulse durations r p as shown in the 
figure, for T = 10 K (a is the rotation angle on the Bloch sphere). Dotted line shows unperturbed 
oscillations. The wavefunction localization widths are l e = 4.9 nm, = 4.0 nm, l z = 1 nm. 
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Figure 5: The final QD occupation after two ir/2 pulses separated by time interval At for pulse 
durations as shown. 



of the optically controlled exciton dynamics for r p = 1 ps (the characteristics for other durations is 
easily obtained by scaling). According to (55-57), the overlap of these spectral characteristics with 
the phonon spectral density gives the perturbation of the coherent carrier dynamics. 

All the spectral densities Ri(co — t{) contain contributions from all the phonon branches and 
all carrier-phonon coupling channels. For the LO phonons, they are peaked at u = ei + £1 (at 
low temperatures). One should note that the denominators £1 and Q + ei are of the same order 
of magnitude, while the strength of the spectral density is much stronger for I > 2 because of the 
charge cancelation effect decreasing the ground-state contribution. Therefore, the response form the 
higher states is stronger. Physically, it corresponds to phonon-assisted excitation of a dark state. 
Obviously, the spectral functions 5 a _ c (u;) must extend to this high frequency sector. 

For acoustical phonons, in contrast, the spectral density are concentrated at low frequencies, so 
that the denominators (cu + ei) 2 for I > 2 are at least two orders of magnitude higher than the typical 
frequency uj 2 appearing in the I = 1 term. Therefore, the contribution from the acoustic phonons 
is restricted to the ground state term. For the piezoelectric coupling, this term is strongly reduced 
by the charge cancelation; therefore, the effect of this coupling on the exciton coherence may be 
neglected as long as the electron and hole densities for the excitonic ground state overlap. 

Damping due to the LA phonons. Let us first discuss the case of relatively long pulses (picosec- 
ond durations) and discuss the contributions from the LA phonons only, neglecting the existence of 
the higher states. The results of such calculations is shown in Fig. 4 for Gaussian pulses [r p is the 
full width at half maximum of the pulse envelope f(t)]. 

The oscillations are almost perfect for very short pulses (~ 1 ps), then loose their quality for 
longer pulse durations (~ 10 ps). Although this might be expected from any simple decoherence 
model, the striking feature is that the effect dramatically grows for higher oscillations, despite the 
fact that the whole process has exactly constant duration. Even more surprising is the improvement 
of the quality of oscillations for long pulses (~ 50 ps) where, in addition, the first oscillation is nearly 
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perfect. 

For growing number of rotations n, the nonlinear pulse spectrum S a (uj) develops a series of max- 
ima of increasing strength (Fig. 3a,b,c). The position of the last and highest maximum corresponds 
approximately to 27m/r p , in accordance with the semiclassical resonance concept. However, spec- 
tral components are also present at all the frequencies 27rn'/T p , n' < n, which is due to the turning 
on/off of the pulse. It is interesting to note that for high n, the low-frequency part of S a (u) does not 
evolve with n (Fig. 3d). It is now clear that there are two ways of minimizing the overlap: either the 
pulse must be so short that all the maxima of S a (uj) are pushed to the right into the exponentially 
vanishing tail of the reservoir spectral density R(uj), or the pulse must be very long, to "squeeze" 
the spectral function near uo = and thus reduce its area. In the latter case, the maxima developing 
with growing number of oscillations will eventually overlap with R{u) destroying the coherence. 

Although it might seem that speeding up the process is the preferred solution, it is clear that 
this works only because no high frequency features are included into the present model. In reality, 
speeding up the dynamics is limited e.g. by the presence of excited states and non-adiabatically 
enhanced LO phonon coupling (see below). Moreover, it turns out that the resulting dynamics 
is actually not fully coherent. It has been shown [31] that superposition of states created by an 
ultra-short ir/2 pulse becomes corrupted, preventing a second pulse (after some delay time At) 
from generating the final state of (n) = 1 with unit efficiency. In order to prove the fully coherent 
character of carrier dynamics it is necessary to demonstrate the stability of the intermediate state in 
a two-pulse experiment. The simulations of such an experiment are shown in Fig. 5. A short pulse 
(r p = 0.1 ps) creates a superposition of bare states (surrounded by non-distorted lattice) which then 
decohere due to dressing processes [29,14]. As a result, the exciton cannot be created by the second 
pulse with unit probability [31]. For a longer pulse (r p = 1 ps), the lattice partly manages to follow 
the evolution of charge distribution during the optical operation and the destructive effect is smaller. 
Finally, if the carrier dynamics is slow compared to the lattice response times (r p ~ 10 ps), the lattice 
distortion follows adiabatically the changes in the charge distribution and the superposition created 
by the first pulse is an eigenstate of the interacting carrier-lattice system, hence does not undergo any 
decoherence and the final effect is the same for any delay time (its quality limited by decoherence 
effects during pulsing). In fact, splitting the 7r-pulse into two corresponds to slowing down the 
carrier dynamics which, in the absence of decoherence during delay time, improves the quality of 
the final state, as seen in Fig 5. 

Effects of LO phonons. The above discussion focused on the relatively long pulse durations 
when the contribution from the LO phonons is negligible. However, an interesting manifestation of 
the resonance effect may be observed for sub-picosecond pulses, when the excitation of longitudinal 
optical (LO) phonons becomes important (Fig. 6). It should be noted that, although the coupling 
between the ground state and the LO modes is strongly reduced by charge cancelation, effects 
involving higher (dark) exciton states may still have considerable impact [30,53]. First, let us note 
that the interaction energy for LO phonons is comparable with their frequencies which results in a 
pronounced reconstruction of the spectrum and appearance of polaron states [11-14,20] which, due to 
strong non-diagonal couplings, mix various excitonic levels. This is manifested in the redistribution 
of the exciton occupation among different states even for arbitrarily long and weak pulses, where the 
dynamical contribution from the LO phonons vanishes (Fig. 6c). It may be easily verified that in 
the limit of slow rotation, when K StC (ei) — > for I > 1, the "static" occupation of the higher states is 
equal to the usual perturbative correction to the state |1). When S a - C (uj) are sufficiently broad, i.e. 
for very short pulses or for large number of rotations, the terms I > 1 contribute also dynamically. 
The effect depends on the relative positions of the narrow LO phonon features and of the induced 
dynamics frequencies and the occupation of the excited states is a non-monotonous function of the 
pulse duration, due to the oscillatory character of the spectral functions S a _ c (u;) [35]. As seen in Fig. 
6, for t p = 0.09 ps, the frequency of induced dynamics overlaps with the phonon frequencies, leading 
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Figure 6: Pulse-intensity-dependent Rabi oscillations for sub-picosecond pulse durations as shown; 
solid: ground state, dashed: M = 1,7V = 0, dash-dotted: M = 2, N = 0, dotted: ideal oscillation), 
where M denotes angular momentum, N. 



to large damping. For longer durations it shifts towards lower frequencies, decreasing the impact on 
the system dynamics. However, the damping is weaker also for shorter durations (r p = 0.03, Fig. 
6a): e.g. for a = 6tt the largest maximum is now to the right of the reservoir response frequencies 
oj + ei, while for lower values of a the magnitude of the spectral functions S^b.c is smaller. 

Conclusion. The above analysis shows that damping of pulse-area-dependent Rabi oscillations 
due to interaction with lattice modes is a fundamental effect of non-Markovian character: it is due 
to a semiclassical resonance between the optically induced confined charge dynamics and the lat- 
tice modes. The destructive effect may be minimized both by speeding up and slowing down the 
dynamics. However, in the former case, the system passes through unstable (decohering) states. 
Moreover, fast operation on a real system induces many undesirable effects: transitions to higher 
states, bi-exciton generation or resonant LO phonon dynamics. On the other hand, for slow oper- 
ation, the number of "good" oscillations is limited. Thus, it is impossible to perform an arbitrary 
number of fully coherent Rabi oscillations on an exciton confined in a quantum dot. 

By increasing the pulse duration as r p ~ a 2 , the phonon effect on the exciton dynamics may be 
kept constant. However, in this case the achievable number of oscillations is strongly restricted by 
the exciton lifetime and other (thermally activated) processes. Eliminating the radiative losses e.g. 
by using stimulated Raman adiabatic passage instead of a simple optical excitation [40] seems to be 
a promising direction from this point of view. 

The model presented above accounts for the decrease of the quality of Rabi oscillations in the 
short duration range observed in the experiment [5]. It predicts, however, that this trend is reversed 
for longer pulse durations. The quantitative value of 96% for the first maximum of the oscillations 
with r p = 1 ps agrees very well with the experimental result [4] although the following extrema are 
much worse in reality than predicted here. This suggests an increased lattice response at higher 
frequencies which may be due to more complicated wavefunction geometry [22] or electric-field 
induced charge separation leading to strong piezoelectric effects [29]. Nevertheless, as far as it may 
be inferred from the experiment [4], the decrease of the oscillation quality seems to saturate after 
one full Rabi rotation, as predicted by the model calculations. 

7 Optimal control over a QD qubit 

Introduction. As we have already pointed out, any fast change of the state of the carrier sub- 
system leads to spontaneous processes of lattice relaxation that affect the coherence of the carrier 
state. In the previous section we have shown that coherent control is recovered if the evolution of 
the carrier subsystem is slow (adiabatic) compared to the typical timescales of the lattice dynamics. 
Thus, the requirement to avoid traces of the carrier dynamics in the outside world favors slow oper- 
ation on the carrier subsystem, contrary to other decoherence processes (of Markovian character), 
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like radiative decay of the exciton or thermally activated processes of phonon-assisted transitions to 
higher states. The latter have the character of an exponential decay and, for short times, contribute 
an error 5 = r g /rd, where r 9 is the gating time and is the time constant of the decay. Means 
to minimize this contribution by speeding up the control operation have been proposed by select- 
ing materials to provide favorable spectrum characteristics [10] or by applying techniques reducing 
unwanted transitions [54,55]. 

In this section we consider the interplay between these two contribution to the error for the 
solid-state qubit implementation using excitonic (charge) states in quantum dots (QDs) [37], with 
computational states defined by the absence (|0)) or presence (|1)) of one exciton in the ground state 
of the dot, operated by resonant coupling to laser light. We show that it leads to a trade-off situation 
with a specific gating time corresponding to the minimum decoherence for a given operation [56]. 



The system and its model. Since it has been shown experimentally [28] that coherence of super- 
positions induced by short pulses is unstable, it seems reasonable to perform operations on dressed 
states, i.e. on the correctly defined quasiparticles of the interacting carrier-phonon system [57]. This 
may be formally achieved by employing the solid-state-theory concept of adiabatic switching on the 
interaction (as done in [36], cf. [58]) to transform the states of the noninteracting system into the 
states of the interacting one. Thus, we assume adiabatic switching on/off of the interaction with 
phonons by appending the appropriate exponent to the original interaction Hamiltonian, 



Hint = e" 



|i><i|£(^4+^V) 



(62) 



where e = + . Fj~ are the deformation potential coupling constants between the ground excitonic 
state and the longitudinal acoustical phonons (20) (this is the only contributing interaction mecha- 
nism for the timescales discussed here). The operator S now becomes 

S(t) = U c (t)e-^\l)(l\ul(t), 

where the free evolution operator is obtained from (27) by truncation to the two lowest states. The 
general formula (47) may now be used with the bare state |Y>o)- The adiabatic procedure assures 
that it is transformed to the dressed state before comparing it to the density matrix p, so that the 
fidelity is defined with respect to the stable, dressed states. 



The irreversible error. The form of the operator Y{ui) = Y\\{u) is obtained in the same way as 
in the previous section. The only difference with respect to (48) is that the oscillating terms now 
vanish in the long-time limit, due to the adiabatic switching on/off described above. For the present 
purpose it is convenient to write the result in the form 

Y(u) = -F( W )(|1X0| - |0)(1| + |0)(0| - |1)(1|) 

+V(-^)(|0)(1|-|1)(0| + |0)(0|-|1)(1|), 

where 

/OO A 
dre iujT f- e^ T \ 
-oo UT 

Since in the quantum information processing applications the initial state of the quantum bit is 
in general not known, it is reasonable to consider the error averaged over all input states. Let us 
introduce the function 

S(u;) = co 2 \(iPo\Y(u J )\^)\l v , 
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where \ipo) is a state orthogonal to |V>o) and the average is taken over the Bloch sphere. By restricting 
(47) to the two-level case, the average error may now be written 

5 = [ ^R{u)S{u>), (63) 

J UJ Z 

where R(lo) = R { ® P \lu). 

The averaging is most conveniently performed by noting that 

Y(u) = -F(u)\+)(-\ + -F*(-u;)\-}(+\, 

UJ UJ 

where |±) = (|0) ± \l))/y/2. Choosing 



i; = cos-|+)+e^sin-|-), ^ = sin -|+) - cos -|-), 



one gets 

2 

F(u)e^ cos 2 - - F*(-u)e-^ sin 2 - 



S(u) 

which, upon averaging over the angles 0, ip on the Bloch sphere, leads to 



s(uj) = ±-(\f(uj)\ 2 + \f(-uj)\ 2 ). 



12 

Let us now consider a Gaussian pulse for performing the quantum gate, 



/(*) = ~J^ e ~~ 2 



2vTT p 



Here r p is the gate duration, while a is the angle determining the gate, e.g. a = | is the Hadamard 
gate, while a = ir is a x (bit flip). The function |F(u;)| 2 that carries all needed information about 
spectral properties of the system's dynamics may be approximately written as 



|F±H| 2 «a 2 e T K^^ P ) (64) 

As may be seen from (63) and (64), for a spectral density R(u) ~ to n the error scales with the 
gate duration as T p _n+1 and T~ n+2 at low and high temperatures, respectively. Therefore, for n > 2 
(typical e.g. for various types of phonon reservoirs, see section 3) the error grows for faster gates. 
Assuming the spectral density of the form R(u) = Rdp^ 3 for low frequencies [in accordance with 
(25) at T = 0], we obtain from (63) and (64) 

5 = -^a 2 Rdpt~ 2 , at T = 

This leading order formula holds for 5 <C 1. Also, if we introduce the upper cut-off, the error will be 
finite even for an infinitely fast gate (see Fig. 4); this is the ultrafast limit discussed in Section 4. 



Trade-off between two types of decoherence. As we have shown, in our model the error 
grows as the speed of gate increases. This could result in obtaining arbitrarily low error by choosing 
suitably low speed of gates. However, if the system is also subject to other types of noise this 
becomes impossible. Indeed, assuming an additional contribution growing with rate 7m, the total 
error per gate is 

o = —j- + 7mt p , 7nM = T^a Rdp, 7m = — , (65) 
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Figure 7: (a) Combined Markovian and non-Mar kovian error for a a = n/2 rotation on a qubit 
implemented as a confined exciton in a InAs/GaAs quantum dot, for T = (solid lines) and T = 10 
K (dashed lines), for two dot sizes (dot height is 20% of its diameter). The Markovian decoherence 
times are inferred from the experimental data [28]. (b) Spectral density of the phonon reservoir 
R{lo) at these two temperatures and the gate profile S(oj) for a = tt/2. 



where r r is the characteristic time of Markovian decoherence (recombination time in the excitonic 
case). As a result, the overall error is unavoidable and optimization is needed. The formulas (65) 
lead to the optimal values of the form (for T = 0) 

3/2o^Rdp\ 1/3 (2 2p y/ 3 

Omin = ^ I ^2 I ' f ° r T P = ^3" ^DP r r J ■ (66) 

For the specific material parameters of GaAs, the optimal gate time and minimal decoherence 
resulting from Eqs. (66) are 

Tp = a 2/3 1.47 ps, <5 min = a 2 / 3 0.0035. 

The exact solution within the proposed model, taking into account the cut-off and anisotropy (flat 
shape) of the dot and allowing finite temperatures, is shown in Fig. 7. The size-dependent cut-off 
is reflected by a shift of the optimal parameters for the two dot sizes: larger dots admit faster gates 
and lead to lower error. 

It should be noted that these optimal times are longer than the limits imposed by level separation 
[37,54,55]. Thus, the non-Markovian reservoir effects (dressing) seem to be the essential limitation 
to the gate speed. 



Error reduction by a double-dot encoding. As we have seen, for long enough pulses the 
fidelity of a one-qubit rotation depends on the low-frequency behavior of the phonon spectral den- 
sity. In this section we will show that for a different encoding of the quantum logical states this 
behavior can be made more favorable, even without changing to a different material. To be specific, 
we consider a qubit encoded by a single exciton in a double-dot system, with the logical values 
corresponding to the exciton location in one or the other dot. The physical reason for the weaker 
lattice impact at low frequencies is that phonons with wavelengths longer than the distance between 
the dots cannot distinguish between the two exciton positions and, therefore, cannot contribute to 
decoherence. 

As in the preceding subsections, we restrict the discussion to the two states, denoted |0) and 
|1), corresponding to the exciton position in one of the two QDs. We assume that it is possible to 
perform rotations in this space, analogous to the resonantly optically driven rotations of the single- 
dot excitonic qubit discussed above. Then, the qubit Hamiltonian Hq is again the 2-dimensional 
restriction of (5) with Ai = 0. The interaction Hamiltonian is 

tfmt = |0)(0| £ (ffftt + if )* 6fc ) + £ ( F W 6 t + F (U. 6fc ) t (67) 
k k 
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where F^' 1 ^ are the coupling constants for the two exciton localizations. Again, we include only the 
DP coupling and suppress the corresponding upper index "(DP)". Let us assume that the dots are 
placed at z = ±D/2. If we neglect the possible differences between the geometry of wavefunctions 
in these dots then, according to the definitions (9) and (10), the coupling constants differ only by a 
phase factor, 



F, 



(0,1) 



D/2 F k 



where F k are the coupling constants for an exciton located at the origin. 
With the help of the Weyl operator 



W = exp 

we define the new bosonic operators 



7fc 



h = Wb k W^ = b k - 7fe . 
In terms of these, the interaction Hamiltonian has the form (up to a constant) 



int 



-e\t\ 



|i><i|£(i^l + W* 



where 



F k 



h(0) 



k z D 
2i sin Ffc. 



(68) 



(69) 



(70) 



(71) 



The physical meaning of the above manipulation is that the lattice excitations are now defined 
with respect to the new equilibrium, corresponding to the presence of the exciton in one of the dots. 
Again, we have introduced the factor describing the adiabatic switching on/off of the carrier-phonon 
interaction. 

The interaction Hamiltonian is now formally identical to that used in the previous subsection, 
but the spectral density now is 

R(u) = i? D p^ (^j u 5 [n B {uj) + l]g(u), 

where g(0) = 1. The non-Markovian error for sufficiently long gate durations at T = may now be 
estimated as 

1 „_ 1 /m 2 

' 'p 



12 3 \ ci , 

Compared to the corresponding value for the simple encoding (65), this shows much faster decrease 
with growing pulse duration. By combining this error with the Markovian decoherence rate resulting 
from the final exciton lifetime and optimizing with respect to the pulse duration one finds 



2 D 



1/5 



for r n 



l a2RBP (| 



1/5 



The comparison of the error with single- and double-dot encoding at T = 10 K is shown in Fig. 8, 
with the exciton decay time the same as in the previous subsection. For the double-dot encoding, 
the error is lower by a factor of 3 and the optimal gating time is shifted to lower values. 
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Figure 8: (a) Combined Markovian and non-Mar kovian error for a a = tt/2 rotation on a qubit 
implemented as a confined exciton in a double InAs/GaAs quantum dot separated by D = 5 nm for 
T = 10 K (solid line), and for the simple encoding (dashed line), (b) The spectral densities of the 
phonon reservoir at T = 10 K (solid: D = 5 nm, dot-dashed: D = 15 nm, dashed: single dot) and 
the nonlinear pulse spectrum (dotted). The wavefunction width was taken to be 4.5 nm in-plane 
and 0.9 nm in the growth direction. The difference between electron and hole localization widths 
has been neglected. 

8 Fidelity of a STIRAP qubit 

Introduction. As we have shown in the preceding sections, the fidelity of quantum information 
processing schemes implemented on orbital (excitonic) degrees of freedom is limited due to the 
finite exciton lifetime, usually of order of 1 ps [28,50], and to the adiabaticity requirement excluding 
unlimited speed-up of the operation. On the other hand, spin-based proposals [59], favored by the 
long lifetime of electron spin [60] also suffer from serious difficulties: the spin switching time in 
typical structures is very long due to weak magnetic coupling. It seems therefore natural to seek for 
a scheme in which the logical values are stored using spin states, while the operations are performed 
via optical coupling to the charge degrees of freedom [38,42,39,41]. 

Recently, it was proposed [40] to encode the qubit into spin states of a single excess electron in 
a QD and perform an arbitrary rotation [61] by employing the stimulated Raman adiabatic passage 
(STIRAP) to a state localized spatially in a neighboring dot [62]. Although this passage requires 
coupling to a charged exciton (X~ , or trion) state which has a finite lifetime, this state is never 
occupied (in the ideal case) so that the scheme is not affected by the decoherence resulting from its 
decay. 

From the discussion presented in the previous sections of this chapter it is clear that the solid state 
QD structures, where the new implementation of these quantum-optical schemes is proposed, differ 
essentially from the atomic systems, where these procedures are successfully applied [63]. It may 
be expected that the coupling to the lattice modes will play an important role in such sophisticated 
quantum control schemes, restricting the range of parameters, where the coherent transfer may be 
performed with high fidelity. 

In these section we briefly discuss the phonon effects on a STIRAP qubit. The complete discus- 
sion exceeds the scope of this chapter and may be find in [64] . 

The Raman adiabatic passage. The arbitrary rotation of the spin qubit between the states |0) 
and |1) (different spin orientations in one dot) may be performed with the help of an auxiliary state 
1 2) [61] (electron in another dot). All these three states are coupled to a fourth state (trion) |3) 
by laser beams Oo>^i>^2- In order to achieve the Raman coupling, the detunings from the corre- 
sponding dipole transition energies must be the same for all the three laser frequencies. Therefore, 
we put uj n = £3 — e n — A, n = 0, 1, 2. The envelopes of the first two pulses are proportional to each 
other, 

7T 

Mo(t) = ^oi(*)cosx, Qi{t) = Q i(t)smx, X G (0, -). 
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Upon transition to the "rotating" basis \n) = e % ^ nt Sn ^\n), 5 n = 5 n — So, n = 0,1,2 and 
suppressing the constant term, the RWA Hamiltonian may be written 

H C = A|3><3| + ^noi(t)(|B><3| + |3>(B|) (72) 



+\tl 2 {t) (e^|2)(3|+e-^|3)(2|). 



2 

where 

\B) = cosx|0) + e iSl sinx|l), \D) = - sm X \0) + e iSl cosx|l), 

Hence, the pulses affect only one linear combination of the qubit states, the coupled (bright) state 
\B), while the other orthogonal combination, \D), remains unaffected. The Hamiltonian (72) has 
the eigenstates 

\a ) = cos6\B) -e^ 2 sin0|2), (73) 
|a_) = cos<P(sme\B) + e l52 cos9\2}) -sin^|3), (74) 
|a+) = sin^sin^S) + e^ 2 cos0|2)) + cos^>|3), (75) 



where 



1/2 



n 1 / A 

tan^ = — , sincp = —= 1 — 

The corresponding eigenvalues are 

A = 0, A± = i ± ^A 2 + ^ 2 + ^ . (76) 

The system evolution is realized by adiabatic change of the pulse amplitudes (in this application, 
the detuning remains constant). Initially (at the time s), both pulses are switched off, hence <f> = 0, 
then 0,2 is switched on first, hence also = 0. Therefore, \ao) coincides with \B) and |a_) with 
1 2). During adiabatic evolution of the parameters, the states move along the corresponding spectral 
branches. As shown in Ref. [61], performing the transfer from 9 = to = ir/2 and than back with 
a different phase 8% of the Q.2 pulse results in a rotation in the qubit space |0), |1) around the axis 
determined by x an d by the relative phase 5\ between an d fii- The rotation angle is equal to 
the difference of the 5% phases in the first and second pulse sequence. Ideally, the state |2) is only 
occupied during gating, while the state |3) is never occupied. 

The above procedure works under assumption that the evolution is perfectly adiabatic. However, 
any change of parameters can never be infinitely slow and the probability of a jump from \oq) to one 
of the two other states \a±) remains finite, leading to non- vanishing occupation of the trion state 
and to decoherence. In order to avoid this error, one has to impose the usual adiabaticity condition 

|A±| » l/r , (77) 

where tq is the duration of the process. 



Phonon induced decoherence. In the presence of phonons, the fundamental adiabaticity con- 
dition (77) is supplemented by the additional requirement to avoid phonon-assisted processes. It 
may be shown [64] that the only phonon coupling term that contributes at the leading order of the 
perturbation is 

H int = \2)(2\J2P k (b k + bl k ), (78) 
fc 
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Figure 9: (a) The contributions to the spectral density at T=0: DP coupling to LA phonons (solid) 
and piezoelectric coupling to TA (dashed) and LA (dash-dotted) phonons. (b) The high-frequency 
behavior of the DP contribution. 

where is obtained in the way analogous to (71) by shifting the phonon modes to the equilib- 
rium appropriate to the occupation of the first dot (states |0),|1)). This time, however, coupling 
constants for a single electron must be used. Vanishing of all the other contributions results from 
the indistinguishability of spin states by phonon interactions, from the large mismatch between the 
trion creation energy and phonon energies and from the fact that the state |3) is not occupied in the 
ideal case. Since now the single, uncompensated charge carrier is shifted between different spatial 
locations, one may expect a considerable contribution from the piezoelectric coupling to acoustical 
phonons. Indeed, as shown in Fig. 9, this coupling dominates at low frequencies, while for high 
frequencies it decreases rather fast due to vanishing geometrical factors (14-16) in the strongest 
confinement direction. In the high-frequency sector, the deformation potential coupling dominates, 
with the oscillatory tail characteristic to a double-dot structure. 

Even though H mt ((78)) is diagonal in the bare carrier states, it still gives rise not only to the 
pure dephasing effect discussed in the section 7 but also to transitions between the trapped carrier- 
field states (73-75). The probability of these phonon-induced transitions becomes very high if the 
spacing between the trapped energy levels falls into the area of high phonon spectral density and 
the contribution to the error resulting from such transitions is approximately proportional to the 
process duration, 

5 ~ R(X±)r , (79) 

with some additional effects appearing due to the pure-dephasing broadening of the A± levels if they 
are placed in the narrow local minima in the tail of R(u). These strong decoherence processes may 
be avoided by either decreasing the trapped level separation (low-frequency regime, exploiting the u> n 
behavior of spectral densities for u — ► 0) or increasing it beyond the cut-off (high frequency regime). 
In both cases one encounters a trade-off situation, due to the opposite requirements for phonon- 
induced jumps (short duration) and for the fundamental adiabaticity condition and pure dephasing 
(slow operation): In the low-frequency regime, avoiding phonon-induced transitions contradicts the 
condition for avoiding non-adiabatic jumps between the trapped states, which may be overcome only 
by considerably extending the process duration. In the high-frequency case, there is competition 
between the pure dephasing and the phonon-induced transitions that is overcome by increasing the 
trapped state splitting, taking advantage of the particular structure of the phonon spectral density 
for a double dot structure. 

In order to provide some quantitative estimations of the error, in Fig. 10 we present the result of 
the numerical optimization of the pulse parameters for |A_| located at the minima of R(u) around 
6 ps -1 and 12 ps -1 [64]. It turns out that a very high fidelity may be achieved for reasonable 
pulse parameters and for acceptable duration of the control sequence (note that tq is a measure of 
duration of a single transfer; taking into account two transfers which should be well separated in 
time leads to the total duration which is roughly ten times longer). 
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Figure 10: The minimal achievable error as a function of the transfer duration corresponding to the 

optimal pulse parameters with both |A_| in the minimum of R(u>) at 6 ps -1 (a) and at 12 ps~ x (b). 

9 Conclusions 

We have discussed phonon-related corrections to externally controlled quantum coherent dynamics 
of carriers confined in quantum dots. We have shown that the carrier-phonon interaction leads to 
lattice re-laxation (dressing process) after a change of the carrier state. This relaxation process leads 
to entanglement between the carrier and lattice subsystems and thus to decoherence of the carrier 
states. 

In order to discuss the possible reduction of this decoherence effect by slowing down the carrier 
dynamics, we have developed a perturbative method for calculating the phonon corrections to ar- 
bitrary evolution. This method allows one to qualitatively predict the phonon effect based on the 
nonlinear spectral properties of the induced carrier dynamics and on the general knowledge of the 
spectral properties of the lattice modes. 

With this approach, we have studied the Rabi oscillations of exciton occupation in a QD. We 
have shown that the lattice response is resonantly driven by a combination of the (linear) pulse 
spectrum and the Rabi frequency. It turns out that the quality of the oscillations improves both for 
fast and for slow dynamics, but only in the later case the intermediate states are stable. 

Applying the perturbative procedure to a QD implementation of a quantum bit we have shown 
that the error for a single-qubit operation may be decreased when one approaches the adiabatic 
limit, in which the lattice modes reversibly follow the carrier evolution. This requires, however, long 
process durations, which increases the contribution from other decoherence mechanism. As a result 
of the interplay between these two contributions, a specific pulse duration appears which optimizes 
the fidelity of the process. The error may be slightly reduced if the simple qubit encoding in a single 
QD is replaced by a double-dot encoding, with qubit values defined by the exciton position in one 
of the dots. 

Finally, we have discussed the phonon impact on the rotation of a spin qubit performed optically 
by adiabatic Raman transition. It turns out that by a suitable choice of parameters, such a procedure 
may be performed with very a high fidelity in a reasonable time. 

The presented results show that phonon-related decoherence processes form the essential limita- 
tion to the coherent manipulation of confined carrier states in quantum dots. However, by employing 
more and more sophisticated control techniques the resulting error may be substantially reduced. 
With the help of theoretical modeling presented in this chapter it is possible to optimize both the 
system properties and the operational parameters in order to maximize the fidelity of quantum 
coherent control. 
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